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Abstract. Density functional theory stems from the Hohenberg-Kohn-Sham-Mermin 
(HKSM) theorem in the grand canonical ensemble (GCE). However, as recent work 
shows, although its extension to the canonical ensemble (CE) is not straightforward, 
work in nanopore systems could certainly benefit from a mesoscopic DFT in the CE. 
The stumbling block is the fixed N constraint which is responsible for the failure in 
proving the interchangeability of density profiles and external potentials as independent 
variables. Here we prove that, if in the CE the correlation functions are stripped off 
of their asymptotic behaviour (which is not in the form of a properly irreducible n- 
body function), the HKSM theorem can be extended to the CE. In proving that, we 
generate a new hierarchy of iV-modified distribution and correlation functions which 
have the same formal structure that the more conventional ones have (but with the 
proper irreducible n-body behaviour) and show that, if they are employed, either a 
modified external field or the density profiles can indistinctly be used as independent 
variables. We also write down the A'^-modified free energy functional and prove that the 
thermodynamic potential is minimized by the equilibrium values of the new hierarchy. 
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1. Introduction 



Density functional theory (DFT) has become, since the seminal work of Hohenberg, 
Kohn and Sham and its extension to the grand canonical ensemble (GCE) by Mermin 
0, one of the more widely used techniques in condensed matter theory. It has been used 
in a great variety of systems [§, ^ as, e.g., uniform and non uniform systems in simple 
P, ^ and general fluids, confined fluids ||, |[|, melting and freezing [|^, interfaces 
I nil, etc. as well as in the calculation of electronic properties of all sort of systems 



This impressive amount of work has been spawned by the Hohenberg-Kohn-Sham- 
Mermin (HKSM) theorem which, roughly speaking, states that, either the external 
potential or the density profile can indistinctly be used as independent variables and that 
the thermodynamic grand potential reaches its minimum value when the equilibrium 
density profile is used. More precisely, it has been proved |13| that any positive density 
can be viewed as an equilibrium density for a system in a suitable external potential. 
Therefore, it is in the foundation of all sort of variational principles in the GCE. On 
the other hand, the failure to implement a straightforward canonical ensemble (CE) 
extension is already well known, referred to in (where a GCE series expansion was 
done and an approximate density profile for the CE obtained), can be traced back to 
the fixed N constraint (see Section 2) and it manifests itself in the failure of having 
an Ornstein-Zernike (OZ) equation. In fact, years ago, Ramshaw explicitly worked 
out the two-body OZ equation for a one component system, showing how the stripping 
of the correlation function h{r) off its asymptotic behaviour solves the problems found 
in a more conventional approach. Another approach to a DFT in the CE, equivalent 
to the one given in |^, is discussed in |15|. In [l^ the existence of an extension of 
the OZ equation for the CE is discussed in the framework of a DFT theory. As the 
GCE is an extension of the CE, one may wonder about the need of an explicit CE 
formulation. The problem addressed by that question is about the equivalence between 
both ensembles that holds true in the thermodynamic limit but breaks down for small 
systems. There is a great body of work done on the properties of fluids confined in very 
small cavities which can shrink down to sizes that hold only a few adsorbed molecules. 
Experiments done in porous glasses (with a mean pore radius r ^ 20 — 30 A) [|T^ have 
shown how such a degree of confinement induces strong changes with respect to bulk 
properties [p!0| , 18 1 as, e.g., critical point shifting 0, capillary condensation [19|, freezing 



point decreasing pO[, etc. Earlier approaches to these highly inhomogeneous systems 
included virial expansions adapted to this situation |22|, calculations of the partition 
function [22, ^ and mean field theory |2|]. Numerical simulation techniques have been 
also used to study adsorption [p5|, phase transitions pOf) P^, etc. Within the 



framework of DFT, the most successful DFT that deals with inhomogeneous systems is 
Rosenfeld's fundamental measure functional (FMF) ^ which, in particular, has been 
used in some of the GCE expansions referred to above ([|l,ll^) in order to obtain a CE 
expansion. It is quite clear, then, the interest of having a mesoscopic DFT in the CE. 



Very recently we have also presented, for one component systems, an extension of 
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the HKSM theorem to the CE. 



In this paper we prove the results put forward in ||2^ and extend them to general 
mixtures. The system we consider is a p components mixture (labeled by greek indices 
and interacting between themselves with two-body potentials), each a component with 
Na particles in a box of volume V and under the influence of one body, species 
dependent, external potentials Va^\x.). We first extend Ramshaw's work \T^] not only 
by considering mixtures but, fundamentally, by proving that by stripping the canonical 
correlation functions off of their asymptotic behaviour an A^-modified set of distribution 
and correlation function hierarchies is generated. Specifically, we prove that they can be 
defined by functional derivatives of the canonical partition function Qn with respect to a 
modified external potential and, moreover, that the ra-body direct and total correlation 
functions are related through n-body OZ equations. The physical idea is the same that 
shows up in almost every textbook on Statistical Mechanics when going from the CE 
to the GCE. There an open subset is considered, here we strip the canonical correlation 
functions off their asymptotic behaviour. The difference is that, for us, this is the 
fundamental step, not an intermediate one. We are then able to prove that this step is 
described, to any n-body order, by functional derivatives with the right formal structure. 
We then prove the HKSM theorem in the CE by introducing an A^-modified free energy 
and thermodynamic potential which is minimized by the -modified density profiles. 
In this way, variational principles formulated in the CE have a firm foundation if the 
A^-modified functions are used in them. Let us emphasize that, as any successful theory 
of liquids needs to consider all the correlation functions (which must have the correct 
/iierarc/iica/ structure), a valid OZ equation for the pair functions h and c is not enough 
to extend DFT to the CE. Therefore, one of the crucial points of this paper is the proof 
of the existence of the A^-modified hierarchies, that they are related by the usual formal 
scheme and are linked by n-body OZ equations. As side remarks, let us mention that i) 
the two-body interaction potential restriction can be also lifted by including more than 
one-body external potentials ^9] but, for the sake of clarity, we will restrict ourselves to 
the two-body interaction potentials case and ii) there exists at least one system (lattice 
gases) where the conventional CE approach fails (due to the fixed topology constraint), 
even in the thermodynamic limit, and it will be discussed in a forthcoming paper. 

The layout of the paper is as follows: in Section 2 we discuss the reasons behind 
the failure of a conventional CE approach, finite size effects in the pair function h{r) 
are discussed in Section 3, while in Section 4 it is proved that, if the (already discussed) 
asymptotic behaviour of the pair canonical correlation functions is removed, then a two- 
body OZ equation in the CE can be formulated. Section 5 discusses both the generation 
of the A^-modified distribution and correlation function hierarchies and of the set of n- 
body OZ equations that relate them while, in Section 6, we introduce an A^-modified free 
energy, write down the variational equations for the density profiles and show that the 
associated A^-modified grand thermodynamic potential is minimized by the equilibrium 
density profiles. This completes the extension of the HKSM theorem to the CE and, 
in this way, variational principles formulated in the CE have a firm foundation if the 
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A^-modified functions are used in them. Section 7 summarizes our results. 



2. Failure of the canonical approach 

The system we consider is a p components general mixture where the species are 
labeled by greek indices a = 1, . . . ,p, each one of them has Na particles (indicated by 
ka = 1, . . . , No) which are under the influence of one-body species-dependent external 
fields l/j^^r) that can be varied at will and the box volume is V. The canonical partition 
function is 

Q({N}, V, T) = ^ J e^"+S^^^. ^''(■'-)d{r} (1) 

Here {rjdenotes the configuration space of all the particles, i.e. d{r} = Ylaka ^i^ka}, 
N! = na=i ^^"^ ^ factor has been absorbed in the hamiltonian and 

fields. In particular, Hq is the system's hamiltonian in the absence of external fields 
(contains all the two-body interaction terms) and $Q,(r) = —(3Va^\r). The one particle 
canonical distribution function is usually defined by [^] 



(1). . .\ /E..^(x-r,Je^o+g--^^.'^-(---)d{r} 
"-a W = ( >^5(x-rfeJ ) = — = (2) 




and can also be written as 



(3) 



Both na\^), the conventional canonical distribution function, and ha\x.), the full 
canonical distribution function (as defined by Lebowitz and Percus [^), are the first 
members of two hierarchies of distribution functions. The s-body functions of these 
hierarchies can be generated by functional derivation 



-«(W) = ^^H? ,Z (X.) (5) 

Here a = (ai,...,^^), {x} = (xi,...,Xs) are the species and coordinate sets and 
the difference between both hierarchies is that in the full one different coordinates can 
coincide, while in the other one can not. The GCE definition is formally the same with 
the grand partition function S replacing Qtv- The truncated correlation functions t and 
t associated to these hierarchies and the direct correlation function c are defined by 
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7(s)/r _ In Q-N _ '^tai..!a,_i(Xi, ■ . ■ ,Xa_i) , . 



riLi '^'^S' (xfc) ^^oi (xs 



p{exc) jg excess free energy. In particular, ta\^) = ia\x.) = na\^) = na\'x.) and 
as 

i£U(xi x^) = n(%,(xi,x,) - <)(xi)<)(x2) = <nxi)nS(x2)/^£U(xi,x,) (10) 

€L(X1 X2) = tgl,(xi X2) + <)(xi)5,,,,5(xi-X2) (11) 

the link with the more usual notation is established. Also, up to now, the independent 
variables are the external fields. As by definition is t^^_^(x, y) = 6na\'x.)/6^\{y) the 
change in the density profiles due to changes in the external potentials is 

Sn^a\^) = E / 4'i(x,y)^<fA(y)dy (12) 

If the density profiles {ria^} and the external potentials {$A}can interchangeably be 
used as independent variables, then an inverse kernel £g"^(z,x) (= 5^r,iz)/6n^a\x.)) 
must necessarily exist such that 

E / f2-^(z,x)tS(x,y)dx = 6,,6{z - y) (13) 

In that case 

5$,(z) = E / tg-^(z,x)5n«(x)dx (14) 

When i^yfa ^ is used to define the direct correlation function c^,a (see equation (|37|) ) then 
the OZ equation follows immediately. Now we can prove that, in the CE, the inverse 
kernel t\^a does not exist, while it does in the GCE. The normalization integrals are 

n(^)(xy)dxdy-/ ^"^^^-^"^^ (16) 

and it is easy to show that 

/n«(y) l + /^S(x,y) dy=<^ 1 (17) 
J I — ^^"i GCE 
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It then follows that 



/^S(x,y)dy 



CE 



5aX + 



dln(n'-a\:x.)/za 
d In zx 



GCE 



Equations (|T3|) and (^) are obviously incompatible in the CE while, in the GCE and 
because of the fluctuations, they are not and the kernel r^l is invertible. In consequence, 
the OZ equation (equations ([37| - |38D ) in terms of the conventional correlation functions 
does not exist in the CE. It is also clear that the existence of an OZ equation is a 
previous condition to the formulation of a DFT but, in fact, the equivalence must be 
established between the existence of the whole set of n-body OZ equations and DFT 
(see Sections 4 and 5). 

3. Finite size effects 

If a A particle is fixed in position y, the conditional probability of finding an a particle 
in X is 

(2)/ N 

n.A(x|y) = ^-y';'^^ =nW(x) l + /^S(x,y) (19) 

This causes a density change that depends on {A, y} which, when x and y are far away 
from each other, can also be written as 

n.A;oo(x I y) = n«(x) + ^!p^Ap,{y) (20) 

opx 

Here Apxly) has the meaning of an unknown proportionality factor which, in fact, 
depends on all the {a, x}; {A, y} variables; it is not the difference in density that arises 
from fixing a A particle in y . Therefore, the asymptotic behaviour of the total correlation 
function can be cast in the form (uncorrelated x-y dependence) 

^Uoo(x,y)= — — ApA(y) (21) 

However, to follow this course of action (deriving with respect to the partial density p\) 
has, at least, two inconvenients: 1) the {a, x} < — > {A, y} exchange symmetry would be 
hidden and 2) the formalism would distinguish between diagonal (hal) and non diagonal 
{h^l^a 7^ A) correlation functions in the sense that the long range behavior of the 
latter cannot be easily separated (in equation (^) the species and field indices would 
be different) and, therefore, the generalization that leads to A^-modified hierarchy of 
distribution functions is far from clear. It can also be mentioned that, while stripping off 
the diagonal components is good enough for obtaining an invertible kernel, we preferred. 



for the reasons mentioned above, not to do that. On the other hand, it is known that 



in the CE, the long range behaviour of the correlation functions is not in the nature 
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of an irreducible two body function, it can be written in a separated variables form. 
Therefore, instead of equations (pO|-pl|) we write (redefining Apxiy)) 



/Sjx,y)=^ii^A,.(y) (22) 
and similarly, for riax-oo- Defining 

1 din Ua^ (x) 



A,, 



r dlnnl^^u) dpc 

J dpa 



j ApA(v)dv 

Both v4a(x) and Bx{y) are normalized to one and /'-aAioolXi y) can be written in two 
possible forms 

^iA;oo(x,y) =f^xA^{x)Bx{y) =T^xB^{x)Ax{y) 
Integrating over x we have rQ,A-BA(y) =rciA^A(y) and so 

um (^^\-v d\nn^^\x) d\nn^x\y) /po^ 

W(^'y)-r«A^^^- (23) 

(1)^ 



dhin\ (u 



■du j ApQ(v)dv 



dpx 

Ta\ depends on {a. A} and on the unknown function Apa(x). An A^-modified canonical 
total correlation function /i^^j (A^-modified as implying that it is different from the 
correlation function h^^\ and that its definition is size-dependent) can then be defined 
by 

^iA(x,y)= (24) 
I Kii^^y) = -1 x = y 

~ (2) 

The A^-modified correlation function /i^^ is the total correlation function stripped off of 
its asymptotic behaviour. Notice that, although this modification seems to affect all the 
h^"^^ values, the excluded volume effects (authentically two-body) are not modified, not 
only for the case x = y but also because, as the stripping is done through a separation 
of variables, this ensures that only the long range behaviour is affected and that the 
irreducible two-body component is not affected by the stripping. Therefore, h)^{^ is the 
correlation function with a truly irreducible two-body behaviour, not h)^x- Lastly, as 
equation (|28|) shows that Vax = Txa, the exchange symmetry {a,x} ^{A,y} is clearly 
preserved. The full truncated correlation function t^^^ associated to h^"^^ is 



4'i(x,y) 



^L'Hx)n«(y)/.i^i 



Aix,yj 



n 



x)(5aA5(x - y) 



iS(x,y)-r 
^S(x,y) 



aX ' 



anW(x)anW(y) 



dpa 



X ^ y 
X = y 



(25) 
(26) 
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In order to write the unkwnown constant T functional of h]^l let us first 



show that dua' {^)/dpa = 1- As Vpa = J na''(x)dx, then dpa/dna\v) 



dn. 



(1) 



[yi)/d Pa = 1. Using this result in equation 



1/V and, so, 



it is obtained that 



tS(x,y)dydx=-r,AV^' 



(27) 



which gets rid of the incompatibility shown in equations ([T3| , [T8D , allows writing the 
unknown constant T^x as 



aX 



1 

172 



5aX + 



nx\y)haxi^^yW 



dx 



f28) 



and to interpret it as a sort of average compressibility. This recovers the classical 



Lebowitz-Percus results in a form suitable for our purposes. We think that it is 
quite possible that the factor F^a provides a formal explanation of the changes observed 
in properties of confined systems with respect to bulk properties [^, |18], |19|, in 
particular when the system is near critical points and thus account for the observed 
deviation in the critical point in confined systems [0 . Moreover, it might be possible that 
the influence of this modification becomes noticeable, even in the thermodynamic limit, 
when the system is near critical points. It is also interesting to study the normalization 
of the A^-modified pair distribution and correlation functions. Writing h 
it is obtained that (equations (0), ([T6| ) and (p!7D) 



'-aA 



"'aX^""'' 



n«(x) I n^)(y)5/.S(x,y)dy = T^xV (29) 

n^a\^) J n^^(y)/^S(x,y)dy = -n«(x)5„A - T^xV (30) 

j nS(x, y)dxdy = Na{Nx - 5„a) - TaxV (31) 
and, using equation (pHj) for Tax, it is concluded that equation (|3TD holds if 

j n«(x)dx j n^^\y)dy = N^Nx (32) 

regardless of the value of / naUxhax- Therefore, the normalization of n^^^ is enough to 
ensure the normalization of nP'^ or, in other words, the norm of the iV-modified pair 
correlation function is undetermined. This free flotation effectively lifts the fixed A^- 
constraint and in the next section it is shown to be the change needed for having a 
two-body OZ equation in the CE. 



4. Two-body OZ equation for the modified kernel 

Now we need to prove that the modified kernel t)^^ (x, y) can be inverted. The 
generalization to higher order functions is done in Sections 5 and 6. and it is needed 
not only for having a well defined A^-modified distribution functions hierarchy, but also 
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for formulating a well defined CE set of OZ equations that operate on arbitrary n-body 
configurations. If the fields are arbitrarily varied, we get, from equations (O) and (E6|) 



Sn^^\^)=Yl J tS(x,y)5<fA(y)dy + 5^r.A^y"nW(y)5$A(y)dy 

A A 

But equation (|) shows that, when only $a varies 

5\nQ = -P5F = j ni')(y)5$A(y)dy 
and, as the partial chemical potential for A particles is nx = V~^dF/dpx)v,T, then 

j n^^\y)5<^x{y)dy = -(3V5^x (33) 



we can write 



(34) 



<5n«(x) = ^ / tS(x,y)5$,(y)dy-5^r„,/3l-5/., 

A A 

and using equation (^) it is obtained that 

5nW(x) =Y,I tti(x,y)5$A(y)dy + 5^ j tti(x, y)/55My 

A A 

= E / ?i(x,y)5$A(y)dy (35) 

A 

where $A(y) is 

^'A(y) = $A(y)+/?/^A (36) 

Therefore, an inverse modified kernel t^,o~^(z, x) exists and complies with equation (^) 
if the modified field $ is used instead of Writing 

^ijl (^)X) = = — c^q(z,x) H - 7^7—— (37) 

oria [y^j (x) 

;;(2)^ 



(which defines the A^-modified direct correlation function Crid) equation (|TB|) becomes 



the OZ equation in the CE where the A^-modified pair correlation functions h)^l and c^^^ 
must be used 

/^S(x, y) - £S(x, y) = E / ^^'nz)5S(x, z)/.S(z, y)dz (38) 

So, up to now, we have seen that the constant N constraint induces a long range 
behavior in the conventional pair correlation function incompatible with being an 
irreducible two-body function and have also shown that, if this behaviour (expressed in 
the separation of variables of equation (p3|)) is dropped off when defining the A^-modified 
pair correlation functions (thus getting a truly irreducible two-body function), then an 
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authentic OZ equation for the A^-modified pair correlation functions is obtained. This 
generahzes Ramshaw's results to pair correlations in a mixture. Let us recall that it 
is essential to show that this scheme gives a coherent view of all the possible n-body 
configurations in a liquid mixture, i.e., it is not limited to pair functions. This is done 
in the next two sections. 



5. The A^-modified canonical hierarchy 



For the sake of clarity we limit ourselves, in this section, to prove that we really have 
another hierarchy of distribution and correlation functions modified in such a way that 
the fixed constraint is properly taken into account and, in the next section the new 
hierarchy is used to prove the extension of the HKSM theorem to the CE. 

We first explicitly write down the modified canonical distribution function h for 
more than two bodies and perform the same sort of analysis as done in Sections 2 and 
3 and then turn to the functional analysis techniques in order to show that, by deriving 
with respect to we can generate higher order distribution functions. In this way we 
not only show that functional analysis techniques are simpler but also obtain, from the 
first method, some additional physical insight. As usual, /?/'^^ will be enough. First of 
all, the definition and normalization equations for the three body functions are 



J2 ~ ''*^«) J2 ~ '^'«) J2 ~ ''"^«) 

^ fca lex f^cx I 

/ J] 5(x - r,j E '^(y - r J E '^(^ - r-j) 

\hoc^lcX lex "^Ti I 

I E ^>-r.jE%-r^)E^(^-r-'^)) 

/ E <^(x-r,J E %-r^)E'^(^-r-j) 



a ^ \ ^ T] 

a ^ \ = Tj 
a — X — T] 



The normalization equations are, in the CE 



J '^alvi^^ Z)dxdydz = Na{Nx - 5ax){Nr, - 5arj - (^Ar,) 

/ nW(z)^il(x,y,z)dz = -(<5., + 5,,)/i2i(x,y) 

j 4^^(y)4^Hz)/iai(^'y'Z)dydz = S^xiSar, + Sxr,) 

J n«(x)4')(y)n«(z)/iS,(x,y,z)dxdydz = NJ^xiSa, + Sx,) 



(39) 
(40) 
(41) 
(42) 
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We can write the A^-modified canonical correlation function h^^^ as 

/^il(x, y, z) = /.i%(x, y, z) + 5/^1% (x, y, z) (43) 

The interesting point to emphasize here is that, as Lebowitz and Percus have 
already discussed, the long range behaviour of an n-body function in the CE ensemble 
depends on how the clustering of the n molecules is. Specifically, and in our case [n = 3, 
long range behaviour), we can have the three molecules very far away from each other 
as well as a 2-molecules cluster and the other one very far of them. It is clear that the 
complexity of the accounting for all the clustering possibilities increases exponentially 
with n and that, in principle, one has to define different constants in equations like (^) 
and (0), one for each kind of clustering. This is more than enough justification to 
follow up the functional differentiation approach instead of the one summarized above. 
However, it is also instructive to show that, if we consider the case of the three molecules, 
every one of them very far from each other, then it is obtained, by using equation ([T8|), 
that 



nW(x)n\^^(yXHz)5C(x,y,z)dxdydz = T^,,V'[1 



V 



V 



(44) 



V ' 1/2 ' 'J/2 ' -5/2 J 

If we can take the thermodynamic limit, the constant Tax^^ is, using equation (^2]) , 
1 

73 



aXrj 



nl^^fx) 



^axi^ari + ^Ar?) 



4^^(y)^r(^)^aA,,(X'y'Z)dydz 



,(1) 



(3) 



dx 



(45) 



This shows that, if the size of the system allows this kind of limit, the long range 
behaviour of n-body configurations is (in the CE) described by (up to) n order 
susceptibilities. Lastly, let us also mention that, working along the same lines that 
led to equation (^), it is concluded that the norm of h^^^ is also undetermined. 

In the functional point of view, the argument is as follows: if the densities and 
fields can be thought of as equivalent sets of independent variables, i.e. given {uq} the 
fields {$a} are a unique functional of them and viceversa, then each one of them can 
be written as a functional Taylor series in the form 



5^x{j)5^r,{,2, 
52$Jx 



-<5$A(y)'5$,(z)dydz+. . . (46) 



A A,jj 



^A ly 

Substituting 5$ in equation 
order, that 



5au,S{y. - u) = ^ / 



and equating powers of 5n, it is obtained, up to second 
5ni'^(x) 5^x{j) 



-dy 



(48) 
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-dydz(49) 



Let us first analyze these equations in the GCE, i.e. both equations (^^Tf ) are 
certainly correct. In that case, equation ( ^Sf ) just leads to the OZ equation. Also, as 
equation (^) holds for all {a, x}, {ti;, u}, {u, v} values and, if we use that (see equations 



])) 



S,,(x,y,z) 



2^(1)/ 



5$A(y)5<l>,(z) 



(50) 



f(3)-l 



x,y,zj 



5n 



(1) 



(y) 



(2)-l, 



X, Z 



)] 



(51) 



(3) / \ 
CaA„(x,y,z) 



in equation ( 
the result is 



e2;'(s,u,v) 



multiply by tia ^(s, x), operate with J dx and use equation 



P(2)- 



-^(s,x)t(^) 



„,^vx, y, z)f2-'(y, u)t(2j-i(z, v)dxdydz (52) 



This is the OZ equation for the three-body configurations and the procedure shows 
how to write OZ equations for n-body configurations, i.e. if the first p-body OZ 
equations are needed, a Taylor expansion of both {n^^} and {$a} is performed, one 
series is substituted into the other and terms up to order p — 1 are equated. These 
equations have been discussed before, see e.g. [Bl, Bl, 13311, and we have followed the 



more straightforward Henderson's method. Notice also that c^^J^ has been generated by 
density differentiation of c^^\. Therefore, exists if and only if all the exist. 

Let us now turn to the CE. If equations (^,^) make sense in the CE, then in equations 
(|^,^) the derivatives must be taken w.r.t. $ and, as a consequence, t is replaced by 
t everywhere and equation (|5^ ) (with t instead of t) is the three-body OZ equation for 
the iV-modified correlation functions. Also, equation (^) shows that, as is customary, 
the A^-modified direct correlation hierarchy is built by deriving with respect to the local 
density. As depends, to all orders, on the lower order derivatives i^''^~^,k < n, 

we conclude that the operators 5/5$ and 6/6n generate the hierarchies of A^- modified 
canonical total and direct correlation functions t and c respectively, with t replaced by t 
and $ by $ everywhere, that the usual formal links between both correlation functions 
are kept and, lastly, let us mention that the missing link of how the direct correlation 
hierarchy starts is discussed in the next section (see equation (|65D) where we identify 
the first member of the A^-modified direct correlation with the density derivative of our 
A^-modified free energy functional. 
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Here we propose an A^-modified free energy and thermodynamic potential functional, 
show that the latter is minimized by the equilibrium density profiles and identify the 
first members of the A^-modified hierarchies completing in that way the proof that we 
really have a new distribution functions hierarchy with its associated truncated and 
direct correlation functions hierarchies and that the HKSM theorem can be extended to 
the CE. 

First of all, let us summarize the standard formulation of DFT in the GCE. It 
starts from the HKSM theorem which states that, in the GCE, not only the external 
potential unequivocally determines the density profile, it is also true that the density 
profile unequivocally determines an external potential |13|. Moreover, the grand 
thermodynamic potential functional is a minimum for the equilibrium density profile. 
This result allows to write it as (following in our notation) 

pn = (3F [{n«}] - /" [$,(x) + n«(x)dx (53) 

F is the free energy functional of the system 

I3F [{n(i) }] = {(3 {Kj, + f/AT - In Pjv)) (54) 

Km is the kinetic energy, Un the interaction energy, P/v the grand canonical probability 
density and the average is a grand canonical one. Let us emphasize that the entropic 
term InP/v depends on all the n-body configurations and, therefore, is responsible for 
the need of a full hierarchy of distribution functions. The variational principle that 
determines the density profile is obtained from 

«1.0 (55, 

Defining the intrinsic chemical potential by 

A-rnK'iixi^^^^i^ (56) 

equation ( |55| ) becomes 

; x] = $ Jx) + = <! Jx) (57) 

where /Iq is the chemical potential and, separating the ideal and excess contributions to 
the free energy F = F^^'^^ — F^^^'^^ 

f^piid) = J2f ^L'^x) [In {X^n'^aH^)) - 1] dx (58) 
/3/if)=ln(A^n(i)(x)) (59) 
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A is the thermal wavelength and the one body direct correlation function is defined by 

^ _p(exc) 



ln(A='n«(x))-<ljx) 



(60) 



which is the usual equation for the density profile 0. The next direct correlation 
function is 



cS(x,y) 



6c, 



6^xS{^-y) ^'l'L'^(x) 

5n«(y) ni^)(x) 5n«(y) 



(61) 



In this way the distribution function hierarchies are generated and a density 
functional theory in the GCE is established. 

Now let us turn to the CE formulation. Equations (0) and (p!8| ) show that a 
straightforward formulation is doomed: the inverse kernel does not exist. In order to 
complete the extension of the HKSM theorem to the CE let us write an A^-modified free 
energy functional, show that there exists a variational principle in the density profile 
and that a whole set of A^-modifed distribution function hierarchies is generated by this 
functional. Our proposal is 



(62) 



where now both P/v and the average are canonical and the A^-modified distribution 
functions n are used. The A^-modified grand potential functional to minimize is 



f3Q [{n«}] =/3^[{n«}] J nL^Hx)$a(x)dx 



(63) 



Let us first show that Cl is minimum when the equilibrium density profiles are used in 
the average of equation (|63|); 



i3n={/3\H-y ]fi^-\nP,,] ) = y P^({x}) 

and, with any other probability distribution Pj^ ({x}) such that 
P; ({x}) d{x} = jp^ ({x}) d{x} 



d{x} 



the Gibbs-Bogoliubov inequality pO| states that 



(3Q' = / P; ({x}) 



d{x} > pn 



(64) 



So, the variational equations are, separating as in the grand canonical case the ideal 
and excess contributions, 



£«(x) 



/5^4jy— = ln(A='n«(x))-<|.Jx) 



(65) 
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This equation for the density profile has the same form as in equation (]60[). The two 
body direct correlation function is 



cIa (x. y) = Zm'Z = ""Tm. - TTTTY^ (66) 



5n«(y) n(^)(x) 5n«(y) 



and it has the same form as equation (|61|). Therefore, it is clearly seen that for the 
A^-modified thermodynamic potential defined in equation ([63| ) a HKSM theorem can 
be formulated and a DFT for the CE is established. Moreover, the thermodynamic 
potential of equation (|63D shows that the iV-modified formulation of the distribution 
functions allows us to work in a CE as if it were a grand canonical one. Also, if the 
system's size forces us to explicitly consider size effects, then equations like ( p6| , p^ ) 
contain a recipe to include size effects. A point worth of mention is that in the first 
sections of this paper we wrote the density profile as ^ and now, in this section, we 
have used fia^. In fact, both are the same {ua^ = = ta^ = • • •) but, as the CE free 

(2) ~ (2) 

energy and the A^-modified free energy also depend on n)^l and n)^l respectively, they are 
not the same functionals and so, the hierarchies generated by functional differentiation 
are different. Therefore, we use ha^ just for the sake of notational consistency. 



7. Conclusions 



Summarizing, we have proved that, if the long range behaviour of the CE distribution 
and correlation functions is, roughly speaking, left aside, the HKSM theorem and the OZ 
equation can be extended to the CE. As the dismissal of the long range behaviour can 
be interpreted as taking an open set within the CE (the textbook procedure to obtain 
the GCE) this seems to be an eminently logical result. We have also generated, first 
by direct construction and then by functional differentiation, a new set of A^-modified 
distribution and correlation functions that obey the same rules and have the same formal 
structure that the more conventional ones have. The point here is that the chemical 
potential modifies the external field and these new functions have the correct n-body 
irreducible behaviour. It can be mentioned that a price to pay for that is to give up the 
normalization equations for the more than one-body correlation functions. In fact, only 
the density profile is normalized to the total number of particles and the other functions 
freely float. The direct construction method also shows that the long range terms 
dropped off depend on n-body susceptibilities. We have also shown that the functional 
differentiation procedure is the amenable one when seeking a general theory but the 
direct construction method shows the influence of the n-body compressibilities and is, 
in this sense, physically enlightening. It might even be thought that its influence could 
be felt, even in macroscopic systems, near the critical points but more specifric work is 
needed. These iV-modified functions nal..as{^i---^s),ial..as{^i---^s),Ca}...as{'^i---^s) are 
defined by 



(67) 
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6' InQ 



piexc) 



^CaiT.^a,_i(xi...X^_l) 



(68) 



(69) 



These equations have the same structure as equations (||),(H) and (|^), while the direct 
correlation hierarchy is initiated by 



Lastly, let us also mention that, although in the introduction we have emphasized 
the possible use of these results on mesoscopic systems, the next paper in this series will 
deal with a system (lattice gas) that, interestingly enough, needs of these results even 
in the thermodynamic limit. 
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